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We calculate the ground-state expectation value of scalar observables in the matrix formulation of 
the random phase approximation (RPA). Our expression, derived using the quasiboson approxima- 
tion, is a straightforward generalization of the RPA correlation energy. We test the reliability of our 
expression by comparing against full diagonalization in Ohcu shell-model spaces. In general the RPA 
values are an improvement over mean-field (Hartree-Fock) results, but are not always consistent 
with shell-model results. We also consider exact symmetries broken in the mean-field state and 
whether or not they are restored in RPA. 
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I. INTRODUCTION 

Mean-field theory is the starting point for most micro- 
scopic models of many-body systems such as nuclei. In 
fact, for global studies of nuclear properties, such as nu- 
clear binding energies, mean-field theory in the form of 
Hartree-Fock (HF) and related approximations are still 
the only viable microscopic approach [jjj; otherwise one 
turns to semiclassical methods such as liquid drop || and 
Thomas- Fermi ||. 

Unfortunately, mean-field theory ignores particle-hole 
correlations and can break fundamental symmetries such 
as rotational and translational invariance. The next logi- 
cal step beyond the static mean-field is known as the ran- 
dom phase approximation, or RPA, which can be derived 
in the small amplitude limit of time-dependent mean-field 
theory, and which implicitly allows for small correlations 
in the ground state ||, ||, [| . 

The main applications of RPA have been to transition 
strengths and excitation spectra. One can also calculate 
the correction to the ground state energy due to RPA 
correlations, which was the basis of a recent proposal to 
use mean-field + RPA to compute global binding energy 
systematics [Q. 

There are other potential quantities of interest, namely 
the ground state expectation value of observables other 
than the Hamiltonian. An important application is the 
root-mean-square radius. Any global fit of nuclear sys- 
tematics involves both binding energies and nuclear radii. 
Indeed, most global fits to binding energies and radii 
with nonrelativistic phenomenological interactions such 
as Skyrme require density-dependent forces, although we 
will not deal specifically with that issue in this paper. 
Further challenges to mean-field calculations are isotopic 
shifts in the charge radius 8, and parity- violating elec- 



tron scattering, which is sensitive to the difference be- 
tween the neutron and proton rms radii [^| . 

In section II of this paper we present an expression for 
the ground-state expectation value of a general observ- 
able in RPA. This expression is derived using the qua- 
siboson approximation, in exact analogy with the RPA 
correlation energy, and is a generalization of the RPA 
one-particle density derived by Rowe [|[ [nj. Because 
the RPA ground state is approximated by a quasiboson 
vacuum, one can read off directly the expectation value 
without having to construct an explicit wavefunction. 

Having an expression is not sufficient. Does it pro- 
duce useful and reliable results? After all, RPA is an 
approximation, and because it violates the Pauli exclu- 
sion principle RPA is not even variational. Despite this, 
rigorous tests of the accuracy of RPA calculations have 
been spotty. For example, although the RPA correlation 
energy has been in the literature for decades, the only 
tests were in toy models Recently we compared 

the RPA correlation energy against exact shell-model re- 
sults |l^] . We found that generally RPA gave very good 
agreement-albeit with some significant failures which re- 
duce the reliability of RPA binding energies. 

We continue this program by computing the ground- 
state expectation value of observables both in exact shell- 
model calculations and in RPA. As discussed below, in 
this paper we limit ourselves to scalar operators. Unfor- 
tunately in Ohuj harmonic oscillator model spaces (R 2 ) 
is trivial, so we instead compute and compare the ex- 
pectation values of a one-body operator: the number of 
particles outside the 0d 5 / 2 orbit; and several two-body 
operators: S 2 , L 2 , Q 2 , and the pairing Hamiltonian 
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The results are found in Section III. In general, RPA 
offers an improvement over the mean-field (HF) expec- 
tation values. This is not always the case, however, and 
even when HF+RPA is closer to the exact shell-model 
results than HF values, the improved agreement is some- 
times still unsatisfactory. (The situation is reminiscent 
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of variational calculations, where one typically gets much 
better agreement in energy than in the wavefunction. 
But again, RPA is not a variational theory.) 

In section IV we turn to another issue of relevance to 
RPA: the restoration of symmetries, such as rotational 
invariance, broken by the mean field. We compute the ex- 
pectation value of J 2 and compare the exact shell-model 
result against the Hartree-Fock and RPA values. If RPA 
"restores" the broken symmetry then one expects that 
either one regains the correct ground state value of J 2 
or gets a much better value than the Hartree-Fock value. 
We find that, although often (J 2 ) is improved by RPA 
corrections, one does not always regain a very good esti- 
mate of the ground state angular momentum. We inter- 
pret this to mean that RPA only approximately restores 
broken symmetries. 

Nonetheless, we consider our results relevant. The 
RPA is not new, but we remind the reader that for global 
approaches to nuclear structure the state of the art is still 
mean-field theory; ab initio methods can only be applied 
to very light nuclei, and shell model diagonalization can 
only be applied piecewise. Consistent calculations using 
HF+RPA throughout the nuclides is conceiveable, how- 
ever, and this paper is part of a larger program to rigor- 
ously test the reliability of RPA as an approximation to 
much larger microscopic calculations. 



II. CORRECTIONS TO A GENERAL 
OPERATOR IN RPA 

Before we derive the expectation value of an operator 
in RPA, we first begin with a brief pedagogical review of 
the derivation of the RPA correlation energy, including 
proper treatment of zero modes that correspond to bro- 
ken symmetries. The RPA expectation value of a general 
operator then follows naturally. 

There are many ways to derive RPA §, §. One ap- 
proach is to expand the energy about the Hartree-Fock 
miminum. Let \HF) be the Hartree-Fock Slater determi- 
nant, and let \V(Z)) = cxp(£ mi Z^c^c^H F) be small 
perturbations from the Hartree-Fock state. We follow the 
usual convention where m, n denote single-particle states 
above the Fermi surface in the Hartree-Fock state, or 
"particle" states, while i,j denote "hole" states below 
the Fermi surface. Then one can expand the energy sur- 
face E(Z) in the vicinity of the Hartree-Fock minimum in 
terms of the particle-hole amplitudes {Z m i} up to second 
order: 



E(Z ) = E H F + z*Az+ ^Bz* + ^z T Bz, 



where 



Anj.mi — (^H F 



c)c 



id C m Ci 



HF^ 



By, 



= (HF [[i?,4c,],4cJ HF 



(1) 

(2) 
(3) 



cause i^HF 



We restrict ourselves to real wavef unctions, so that A, 
B are also real. There is no linear term in (|l|) bc- 

HF^j — by the definition of 

the Hartree-Fock state. This quadratic surface can be 
mapped to a harmonic oscillator, replacing the Z mi by 
boson operators b mi 

Hb = E HF + &Ab+ ~#B# + -bBb, (4) 



which has analogous commutation relations, e.g., 

) , etc. The boson Hamil- 



b nj ,[H B ,V mi ] 



tonian can be written in matrix form, which induces an 
additional constant term: 



H B = E HF -±Tr A +^($ b) 



A B 
B A 



(5) 



We want to put ([|) into diagonal form, the first step of 
which is to solve the well-known RPA matrix equation: 



A B 

-B -A 



= n 



(6) 



The solutions have normalization X 2 — Y 2 = 1. 

Before going on one must treat carefully the case of 
zero modes, fi„ = 0. Zero modes arise in RPA when an 
exact symmetry is broken. Suppose the HF state is de- 
formed. The HF energy will not change as the orientation 
is rotated; this invariance is reflected in E(Z) and resur- 
faces as RPA zero modes. (This is also a good check 
of one's RPA codes. One expects three zero-frequency 
modes for triaxially deformed states, and two zero modes 
for an axisymmetric deformed state: rotation about an 
axis of symmetry takes one to the same, indistiguishable 
state.) One cannot define properly normalized X, Y for 
zero modes. Instead one introduces collective coordinates 
Q v and conjugate momenta .PV@, |^, which satisfy 



ap„ - bp; = 
AQ„ - bq: 



iMMlQv, 
i 



My 



-p. 



(7) 



Here M v is a constant, interpretable as mass or moment 
of inertia but whose value depends on the normalization 
of P, Q, which are only constrained by 



Ql-Pu-Qx-P; = i8 X v 



(8) 



Now one obtains a generalized Bogoliubov transforma- 
tion of the boson operators b, w , not only to quasibosons 
[3, fi\ but also to the zero-mode collective coordinate 
operator Q v and conjugate momentum operator V v : 



A,n A >o 



(9) 
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Application of (^) and its Hermitian conjugate to (||) 
puts it in diagonal form: 



h b = e rpa + J2^xpIpx+ 



V 2 
2M„ 



(10) 



The ground state is the quasiboson vacuum, and one can 
easily read off the RPA ground state energy, which is just 
the zero-point energy: 

E RPA = E HF - ^Tr(A) + -Tr(O). (11) 

One can expand 

A = XOX t +PM" 1 P t , (12) 

where the second term is restricted to zero modes, and 
then obtain an expression equivalent to (O), 



n„=o 



2 

v\ > 



(13) 



an expression which explicitly segregates out the contri- 
bution from zero modes. 

The above derivation of the RPA correlation energy 
-±Tr(A) + |Tr(n) is well-known. We derive the RPA 
expectation value of any generator operator O in exactly 
analogous fashion. Define A, B as in Eqs. (|J) and (j^), 
respectively, by replacing the Hamiltonian by the observ- 
able O, that is, 



iiP) 



HF 
HF 



CjC n , [O, c^Cj] 
[O, c^Cj], c m £i 



hf}, 

HF^ 



(14) 
(15) 



and define o m i = {HFUc^Ci, 0]\HF), which does not in 
general vanish. (Because we restrict the Hartree-Fock 
state to real Slater determinants, A,B and o are real.) 
Please note that in computing A, B, and o one still uses 
the Hartree-Fock state from minimizing the Hamiltonian 
H, not from minimizing O; we are finding the expecta- 
tion value (O) in the vicinity of the minimum of (H). 

The form of the operator in the harmonic oscillator 
bosons is 



Ob = Ohf — — Tr A + ( o, o) 



(16) 



By transforming to the collective quasibosons (|9j), we 
again can trivially read off the quasiboson vacuum expec- 
tation value, which is the RPA ground-state expectation 
value, without explicit construction of a wavefunction: 



Orpa — Ohf 



where 



-Tr A+ -Tr0, 
2 2 



= X T AX + Y AY 



X T BY 



Y T BX. 



(17) 



(18) 



Substitution of Eq. (g) with A, B derived from the 
Hamiltonian immediately regains the RPA binding en- 
ergy (pd]). It is important to emphasize again that the 
X, Y used here are those calculated from Eq. using 
the original A, B matrices (from the Hamiltonian); one 
does not compute X, Y using A, B. 

As before we can rewrite ( |T7j ) into an expression with 
explicit segregation of the zero modes: 



Orpa = Ohf + E J2 ( A 

\((l\>0) mi,nj 

~ 2* ^ t ^ t A m i,nj(P m i,i_iQnj,[j, ~ Qmi.^Pnj,^)- 
Ai(Sl M =0) mi.nj 



(19) 



We have confirmed numerically that (|19[) yields the same 
values as (P~7|) . 

As we will see in the next section, Eq. ([l7]) reduces to 
that derived by Rowe for one-particle densities in spher- 
ical nuclei [g (To). 

If O is a scalar, so that it is invariant under rotation, 
translation, etc., then the above is sufficient. Nonscalar 
observables require a little more thought. In particular, 
consider observables that are angular momentum tensors, 
O k „, with a nonzero rank k. Here /i is the magnetic 
quantum number. If one applies the generalized Bogoli- 



ubov transformation (|9|) then the linear terms in o m j do 
not vanish. The linear terms proportional to the quasi- 
boson operators /3, j3' are transitions to excited states; 
but, for nonscalars, there will be terms also proportional 
to <2„ and to V v . Furthermore, there will be quadratic 
terms, such as ffi (3, (3/3, etc., which manifestly do not 
contribute to the ground state, but also terms quadratic 
in Q, V . These terms arise from rotations (and in the 
more general case, translations, etc.): because they corre- 
spond to zero- frequency modes, they connect only to the 
ground state, but in a different orientation. If O k „ is a 
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spherical tensor, we know how it transforms under a ro- 
tation, and the linear and quadratic terms in Q, etc., are 
simply the linear and quadratic terms for small rotation 
angles. Nonetheless, the issue of nonscalar observables is 
not trivial, and we leave it to future work. 

To summarize a first result of the present work, we 
have extended the usual expansion of the Hamiltonian 
operator in the neighborhood of the HF solution to any 
operator, with the purpose to incorporate RPA corre- 
lations in computation of expectation values of relevant 
observables. In the next section we test the reliability 
of this formula in a non-trivial model, that is, the shell- 
model restricted space using realistic Hamiltonians. 

III. RPA VS. SHELL-MODEL FOR SCALAR 
OBSERVABLES 

In this section we test the RPA expectation value 
found in Eqs. (p"7|), (|l9|), against an exact numerical solu- 
tion computed by diagonalizing a Hamiltonian in a shell- 
model basis. We work in complete Ohu spaces where 
the valence particles are restricted to the single-particle 
states of a single major harmonic oscillator shell, e.g., 
ls 1/ 2-0d 3/2 -0d 5 /2 for sd or IP1/2-IP3/2-O/5/2-O/7/2 for pf 
shell. The cores are inert, 16 O for the sd and 40 Ca for 
the pf shell. In the the valence space we use realistic 
interactions: Wildenthal "USD" in the sd shell Q and 
the monopole-modified Kuo-Brown "KB3" in the pf shell 

As mentioned in the introduction, an operator one 
would like to measure is R 2 . Unfortunately in a OHlu 
shell-model space (R 2 ) is trivial. The operator R 2 has 
two pieces, a one-body piece r 2 , which has a constant 
expectation value in a single major harmonic oscillator 
shell, and a two-body piece r(l)-r(2); but because r is an 
odd-parity operator, the two-body piece is non-zero only 
across major shells. Therefore in this paper we consider 
other one- and two-body operators. 

A. One-body operators 

Using the quasiboson approximation, Rowe derived the 
RPA one-particle densities for spherical nuclei Q [To|| ; 
related formulas have been used to compute the isotope 
shift in calcium ignoring the two-body contribution. 
Let pm — (c\,{Cm) be the number of particles excited 
above the Fermi surface into the particle state M. (M 
is not a magnetic quantum number.) Here A rn i^ n j = 
SijS mM S n M and from (|l|) 

A i 

which is Rowe's result for spherical Hartree-Fock states. 
If the HF state is deformed, however, one has to take 
into account zero modes. The second term in (|19|) can 



TABLE I: 


The number of particles 


in the d 3 / 2 + s 


1/2 orbits 


C = /r/2 


orbit for 48 Ca). * denotes nuclides with 


spherical 


Hartree-Fock states. 






Nucleus 


HF 


RPA 


SM 


20 Ne 


1.60 


1 7^ 


1.60 


22 Ne 


1.49 


1 (\A 


1.47 


24 Me; 


2.13 


1 


2.31 


28 Si 


3.87 


3.77 


3.15 


20 


0.24 


n ap, 


0.54 


21q 


0.06 


0.53 


0.38 


22q* 


0.00 


0.67 


0.54 


24 Q , 


2.00 


2.50 


2.25 


48 Ca*t 


8.00 


7.73 


7.78 


22 Na 


1.51 


1.69 


1.52 


26 Al 


2.29 


2.51 


2.18 


19 F 


1.02 


1.10 


1.52 


21 p 


0.84 


0.56 


0.87 


25 Mg 


2.21 


3.09 


1.96 



be simplified if one sums over all particle states M, that 
is, the total number of particles excited above the Fermi 
surface; in that case one can use Eq. (||) and find 

E p m = E E i*w 1 2 - \ N ™° ( 21 ) 

M A Mi 

where N zero is the number of zero modes, iV zero = 2 for an 
axisymmetric Hartree-Fock state, such as for 20 Ne, and 
=3 for a triaxial state, found for 24 Mg. We find often 
leads to negative occupation numbers! Note, however, 
that ^2 Pm is not an angular momentum scalar, nor even 
a spherical tensor of fixed rank, for deformed states. 

Instead we considered the only scalar one-body opera- 
tors in our system, the occupation numbers of a j-shell. 
Table I tabulates the Hartree-Fock, RPA, and exact shell- 
model (SM) values of n(d 3 / 2 ) + n(si/ 2 )- The results are 
rather poor, except for three cases with spherical Hartree- 
Fock states, 22 0, 24 0, and 48 Ca (for which we instead 
tabulate n(/V/2)). 21 is weakly deformed and also shows 
a reasonable, albeit imperfect, improvement as one goes 
from the Hartree-Fock value to the RPA value. 

To explore this issue further, we took 28 Si and lowered 
the d$/2 single-particle energy; eventually the Hartree- 
Fock state changes from deformed to spherical. The re- 
sults are plotted in Figure |l|. Again we see reasonable 
agreement for the spherical region, but poor agreement 
in the deformed regime. 

To summarize our results for one-body operators: we 
regain, for spherical Hartree-Fock states, Rowe's one- 
particle occupation numbers and get improved values 
over the Hartree-Fock occupation numbers. For de- 
formed nuclides, however, the RPA value is generally 
worse than the HF value. The fault does not appear 
to lie in the corrections due to zero modes; in the next 
section, we will find that the RPA expectation value of 
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FIG. 1: n(d3/2) + nisx/%) in 28 Si as the d$/2 single-particle 
energy is lowered relative to the other single-particle energies. 
The solid line is the (exact) shell-model value, the dotted line 
the HF value, and the dashed line the RPA value. 



J 2 is more accurate in the deformed regime than in the 
spherical regime. 



B. Two-body operators 

We now turn to two-body operators, or more properly 
operators with both one- and two-body pieces. (We in- 
vestigated the pure two-body pieces but found no qual- 
itative differences; the pure two-body pieces performed 
neither better nor worse on the whole than the one-body 
pieces, which here are linear combinations of number op- 
erators.) 

In Table || we show results for S 2 (total spin) and L 2 
(total orbital angular momentum). The RPA expecta- 
tion value is generally a significant improvement over the 
Hartree-Fock value, relative to the exact result. On the 
other hand, the RPA values, while closer to the mark, are 
not in very good agreement with the exact shell-model 
values, and sometimes overcorrect to negative, nonphys- 
ical expectation values (this can happen because RPA 
does not respect the Pauli exclusion principle). 



In Table III we consider the expectation value of the 
pairing interaction, P+P where P* = Ylj, m >o a ],m a i,-m> 
and of Q 2 . We also show the ratio of correlation energies 
fl2f which is a measure of how well the RPA binding en- 
ergy tracks the exact binding energy. There appears to 
be no correspondence: a good RPA value for the binding 
energy does not correspond to a good RPA expectation 
value. In particular, note the single-species (oxygen) re- 
sults, where the RPA binding energy is particularly bad; 
yet for these nuclides (P^P) and (Q 2 ) are very good. 

Again we look at the transition from deformed to 
spherical in 28 Si for (Q 2 ) and (ptp), in Fig. |, which 
clearly shows the RPA values are in better agreement in 
the spherical regime than in the deformed regime. (As it 



TABLE II: Ground-state expectation values of S 2 and I? for 
several nuclei in sd and pf shells. For each observable we 
show the SM, HF and RPA estimates. The nuclides have been 
grouped into even-even, single-species, odd-odd, and odd-A. 



Nucleus 




C'2 

s 






r 2 

L 






HF 


RPA 


SM 


HF 


RPA 


SM 


Ne 


0.35 


0.33 


0.26 


15.90 


-0.25 


0.26 


22 Ne 


1.48 


0.48 


0.88 


16.76 


0.31 


0.88 


24 TV T 

Mg 


1.39 


1.38 


1.03 


20.65 


-1.17 


1.03 


Mg 


2.04 


1.14 


1.45 


18.94 


-0.34 


1.45 


28q; 

Si 


1.62 


1.28 


1.45 


21.50 


-0.75 


1.45 


44 rri • 
ll 


1.03 


0.75 


0.64 


30.34 


-2.48 


0.64 


ll 


2.24 


1.20 


1.36 


29.72 


-2.94 


1.36 


Or 




n nit 

u.yy 


1. lU 


on 77 
zy. / / 


o.oo 


i 7n 
1. IU 


(J 


1.50 


0.45 


0.75 


6.80 


0.92 


0.75 


22^, 


2.40 


-0.15 


1.26 


2.40 


6.36 


1.26 


24 o 

u 


2.40 


-0.27 


1.29 


2.39 


6.06 


1.29 


20 p 


2.00 


1.44 


1.74 


14.21 


8.09 


3.55 


22 Na 


2.20 


1.90 


2.14 


21.32 


9.08 


8.07 


26 Al 


3.14 


1.96 


1.45 


29.56 


20.14 


1.45 


46 y 


2.51 


1.50 


1.36 


35.39 


16.33 


1.36 


19p 


1.09 


0.80 


0.87 


12.61 


4.39 


0.22 


21p 


2.11 


0.76 


1.52 


13.31 


5.60 


6.41 


21 Ne 


1.11 


0.44 


1.00 


17.55 


10.11 


3.22 


23 Na 


2.02 


0.88 


1.15 


18.81 


7.46 


3.93 


25 Mg 


2.04 


0.38 


1.73 


22.56 


11.77 


7.68 



TABLE III: Ground-state expectation values of P'P (pairing) 
and QQ. The final column is the ratio of the RPA correlation 
energy to the shell- model correlation energy, and =1 when the 
RPA binding energy is equal to the exact binding energy. 



Nucleus 




Pairing 






Q-Q 




E H f — Er.pa 

E hf —Esm 




HF 


RPA 


SM 


HF 


RPA 


SM 




20 Ne 


2.99 


5.47 


6.81 


715 


825 


793 


0.75 


22 Ne 


3.99 


7.25 


9.31 


876 


1007 


944 


0.97 


24 Mg 


5.99 


10.14 


11.72 


1167 


1263 


1268 


0.92 


26 Mg 


6.99 


11.51 


14.56 


1001 


1104 


1048 


0.94 


28 Si 


8.99 


12.73 


15.16 


1304 


1389 


1214 


0.90 


20 


2.00 


5.18 


7.25 


257 


353 


339 


1.09 


22q 


3.00 


5.83 


6.20 


163 


277 


270 


1.67 


24q 


4.00 


6.52 


6.58 


122 


194 


191 


1.83 



happens, of the nuclides we investigated 28 Si, while con- 
venient for comparing spherical vs deformed regimes, is 
the only nuclide for which the RPA value of Q 2 is worse 
than the HF value, using the original Wildcnthal single- 
particle energies.) This is not universal behavior; as seen 
m Table |T| and will be seen in the next section for J , 
the RPA expectation value for some operators is better 
in the deformed regime. 
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FIG. 2: (Q 2 ), {P f P} in 28 Si as the d 5/2 single-particle energy 
is lowered relative to the other single-particle energies. The 
solid line is the (exact) shell-model value, the dotted line the 
HF value, and the dashed line the RPA value. 



IV. RESTORATION OF BROKEN 
SYMMETRIES? 



TABLE IV: Ground-state expectation value (J 2 ) for selected 
nuclides, grouped into even-even, single-species, odd-odd, and 



odd-A. 


nJllf*lpi1Q 


HF 


RPA 


SM 


20 Ne 


1 (S OR 


-fl 4^ 


n 

u 


22 Np 


17 17 


1 1 fi 

- J — LU 


n 


24 Me- 
lvi 5 


20.13 


-2.52 


o 




18.61 


-1.72 





28 Si 


20.89 


-1.99 


o 


44 rpj 


31.65 


-3.10 


o 


46 rj^ 


31.53 


-5.00 


o 


48 Cr 


29.37 


4.72 


o 


20 q 


6.07 


1.76 





22 q 


0.00 


7.99 





24 q 


0.00 


7.38 


o 


20p 


18.46 


12.41 


6 


22 Mo 

in a 


op; p.7 


i a p;7 
14.0 / 


1Z 


26 Al 


35.98 


27.92 





46 y 


39.56 


20.00 





19 F 


15.12 


5.52 


0.75 


21p 


15.51 


9.47 


8.75 


21 Ne 


19.05 


12.68 


3.75 


23 Na 


19.42 


11.87 


3.75 


25 Mg 


23.87 


14.51 


8.75 



The random phase approximation respects broken 
symmetries by separating out exactly, as zero modes, 
spurious motion. This is sometimes interpreted as an 
"approximate restoration of the symmetry" The 
restoration cannot be exact, because the RPA wavefunc- 
tion is valid only in the vicinity of the Hartree-Fock state 
[Jk| and cannot be extrapolated to, for example, large 
rotation angles. 

Still, we now have a tool to further explore symmetry 
restoration, by computing Casimir operators of symme- 
try groups. Specificially, we consider (J 2 ). Ideally, if the 
RPA restores a broken symmetry, one might imagine that 
one either regains the exact ground state value of (J 2 ) or 
gets very close to it. 



We present our results in Table IV . The pattern is the 
same as with other operators: ( J 2 )is generally better in 
RPA than in Hartree-Fock but not always very close to 
the exact shell-model value. Even worse are the cases 
with a closed shell in HF, such as 22,24 0: the HF value is 
correct, while the RPA value is terrible! 

To examine this issue more closely, in Fig. ^| we again 
plot, for 28 Si, (J 2 ) versus the ^5/2 single-particle energy 
through the transition from deformed to spherical HF 
state. The results are better for the deformed HF state, 
although we obtain slightly negative, and thus nonphys- 
ical, values of (J 2 ). 

An additional test of symmetry restoration would be 
computation of the expectation value of a nonscalar ob- 
servable, such as the magnetic dipole moment or electric 




"(15/2 



(MeV) 



FIG. 3: (J 2 ) in 28 Si as the d 5 / 2 single-particle energy is low- 
ered relative to the other single-particle energies. The solid 
line is the (exact) shell-model value, the dotted line the HF 
value, and the dashed line the RPA value. 



quadrupole moment, for a deformed nucleus with a J = 
shell-model ground state. We have preliminary, unpub- 
lished calculations which suggest that indeed the RPA 
ground state of even-even nuclides retains a significant 
quadrupole moment, another piece of evidence that sym- 
metry is incompletely restored. 

There are other observables one would like to compute 
relevant to broken symmetries. Aside from the Casimir 
itself, the dispersion of a Casimir would be a useful mea- 
sure. For example, consider quasi-particle RPA (QRPA), 
where particle number is broken in the Hartree-Fock- 
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Bogoliubov state. One would like to see the QRPA value 
of the dispersion (N 2 ) — (A) 2 move close to zero. An- 
other example would be proton-neutron RPA (pnRPA) 
or pnQRPA, allowing protons and neutrons to mix, so 
that T z = Z — N is no longer an good quantum number; 
this is applicable to beta decay. Here one might consider 
(T 2 ) — (T 2 ) 2 . Unfortunately, we suspect that the disper- 
sion would also signal incomplete symmetry restoration. 

V. SUMMARY AND CONCLUSION 

We derived a expression for the ground-state expec- 
tation value of observables in the matrix formulation of 
RPA, and tested it against exact shell-model calculations 
for selected scalar operators. The RPA value was in gen- 
eral an improvement over the Hartree-Fock value, but 
failed to be a consistent and reliable estimate of the ex- 
act expectation value. Nonetheless this work should be 
considered a starting point for any modified RPA calcu- 



lations, such as renormalized RPA, etc. 

In particular we considered the expectation value of J 2 . 
If one starts with a deformed Hartree-Fock state, which 
breaks rotational invariance, the RPA approximately re- 
stores rotational symmetry, as evinced by better values 
of (J 2 ). The results are not wholly satisfactory, however, 
as (J 2 ) can take on unphysical (negative) values; fur- 
thermore, if one starts from a Hartree-Fock state with 
good symmetry, the HF value of (J 2 ) is correct while the 
RPA value is large and positive, a disappointing result. 
Thus, while the RPA respects or identifies broken sym- 
metries exactly, one can only characterize the restoration 
of symmetry in the RPA as approximate and somewhat 
unreliable. 
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